A two-mode squeezed light based on a double-pump phase-matching geometry
He Xuan-Jian1, Jia Jun1, Jiao Gao-Feng1, Chen Li-Qing1, †, Yuan Chun-Hua1, 3, ‡, Zhang Wei-Ping2, 3
State Key Laboratory of Precision Spectroscopy, Quantum Institute for Light and Atoms, Department of Physics, East China Normal University, Shanghai 200062, China
School of Physics and Astronomy, and Tsung-Dao Lee Institute, Shanghai Jiao Tong University, Shanghai 200240, China
Collaborative Innovation Center of Extreme Optics, Shanxi University, Taiyuan 030006, China

 

† Corresponding author. E-mail: lqchen@phy.ecnu.edu.cn chyuan@phy.ecnu.edu.cn

Project supported by the National Natural Science Foundation of China (Grant Nos. 11974111, 11474095, 11874152, 11604069, 91536114, 11654005, and 11234003), the Fundamental Research Funds for the Central Universities, China, the Science Foundation of Shanghai, China (Grant No. 17ZR1442800), and the National Key Research and Development Program of China (Grant No. 2016YFA0302001).

Abstract

We theoretically investigate the frequency-nondegenerate and frequency degenerate squeezed lights with a four-wave mixing process (4WM) driven by two pump fields crossing at a small angle. Different from a 4WM process driven by a single pump field, the refractive index of the corresponding probe field, np, can be converted to a value that is greater than 1 or less than 1 by an angle adjustment. In the new region with np < 1, the bandwidth of the gain is relatively large due to the slow change in the refractive index with the two-photon detuning. In this region with an exchange of the roles of the pump and probe beams, the frequency degenerate and spatial nondegenerate twin beams can be generated, which has potential application in quantum information and quantum metrology.

1. Introduction

Quantum enhanced metrology is the use of quantum techniques to improve measurement precision, which has been received a great deal of attention in recent years.[13] The Mach–Zehnder interferometer (MZI) and its variants can provide high precise measurements, based on which the gravitational waves were observed.[4,5] The sensitivity of these interferometers is limited by the vacuum quantum noise injected from the unused port. To further improve measurement sensitivity, Yurke et al.[6] introduced a new type of interferometer in which two four-wave mixers or parameter amplifiers (PAs) occupy the place of two linear beam splitters (BSs) in the traditional MZI. It is also called the SU(1,1) interferometer because it is described by the SU(1,1) group, as opposed to SU(2) for BSs. Because it can be used to improve measurement accuracy, these types of interferometers have received extensive attention both experimentally[716] and theoretically.[1730] Kong et al.[18] firstly proposed a simplified scheme of PA + BS, and it can also beat the standard quantum limit of phase-measurement sensitivity by a similar amount for the SU(1,1) interferometer. The improvement comes from a result of noise reduction of the interferometer because the quantum noise of the outputs of the parametric amplifier is entangled and destructive interference at the beam splitter leads to the quantum noise cancellation. The scheme of PA + BS can reduce the difficulty of the experimental implementation. However, the scheme of PA + BS requires the two outputs of the PA to be frequency degenerate and spatial modes separation. Therefore, it is worth initiating a study on how to generate a two-mode squeezed state of frequency degenerate and spatial modes separation.

The first experimental demonstration of squeezed states of light by Slusher et al.[31] was based on four-wave mixing (4WM) in sodium vapor. Since then, many techniques for producing different types of squeezing have been explored, each with their own advantages and limitations for particular applications.[32] Nondegenerate 4WM in a double-Λ scheme[33] was identified as a possible scheme to generate a squeezed state or squeezed twin beams, as described in Refs. [3441].

The generated twin beams by the 4WM process in an atomic system with higher squeezing degree were firstly realized by McCormick et al.[40,41] based on degenerate pump fields. A single linearly polarized pump beam is crossed at a small angle with an orthogonally polarized much weaker probe beam. The 4WM process amplifies the probe and generates a quantum-correlated conjugate beam, on the other side of the pump (at a higher frequency). In this case, a pair of photons of the (single) pump is transformed, via the 4WM process, into a photon in the probe beam and a photon in the conjugate beam. By modulating the involved ground (excited) state with one (two) laser beam (beams), the gain and squeezing degree can be enhanced.[42,43] The best initial results for two-mode intensity-difference squeezing at low frequencies seem to be ≈ 1.5 kHz[44] to the recently reported ≈ 700 Hz[45] or even ≈ 10 Hz.[46] The generated entanglement between the probe and conjugate beams can realize quantum imaging.[47,48] The cascaded 4WM can generate the quantum correlated triple beams[49,50] and can also be used to realized SU(1,1) interferometers for highly sensitive phase measurements.[16,51] This 4WM process also supports many spatial modes, making it possible to amplify complex two-dimensional spatial patterns.[5255]

Recently, a new 4WM process driven by two pump fields at a small angle was realized,[56] where we only demonstrated the phenomena and results, without giving detailed physical explanations. Turnbull et al.[57] developed a good theory to illustrate a 4WM process driven by a single pump field. In this paper, we describe the two pump field phase matching that can be established between the np > 1 and np < 1 by the angle adjustment, which is different from the single pump case, and describe the large bandwidth of the gain varying with the two-photon detuning. The theoretical range of phase-matching angles for achieving different regions is given.

Our paper is organized as follows. In Section 2, we describe the generation of two-mode squeezed light with two pump fields crossing at a small angle. In Section 3, the phase-matching angle range is given under the condition of two pump fields driven. In Section 4, the gain of two-mode squeezed light driven by one pump field and two degenerate pump fields is compared, and by an angle adjustment the phase matching between the np > 1 and np < 1 is described. Finally, we conclude with a summary of our results.

2. Theoretical model for a two-mode squeezed light

The generated frequency non-degenerate squeezed light by the 4WM process is shown in Fig. 1(a). Based on the degenerate pump fields, the twin beams with higher squeezing degree were generated,[40,41] where the phase-matching geometry is shown in Fig. 2(a). In our scheme, one unique feature of this 4WM process is that two pump fields are applied at a small and adjustable angle, where the phase-matching geometry is shown in Fig. 2(b). The approach to generate the frequency degenerate twin beams is based on the idea of inverting the configuration from Fig. 1(a) to Fig. 1(b). For the geometry of a single pump field shown in Fig. 2(a), with a simple exchange of the roles of the pump and probe beams, single-mode squeezing is generated.[58] To obtain the twin beams with the same frequency that are spatially nondegenerate, the geometry of two non-collinear pump fields shown in Fig. 2(b) is needed.

Fig. 1. Double-Λ scheme on the D1 line of 85Rb configuration for the generation of (a) the frequency non-degenerate and (b) frequency degenerate squeezed lights. The state |g,m⟩ involves the hyperfine levels |5S1/2, F = 2,3⟩. The hyperfine splitting of the excited state is not resolved due to Doppler broadening.
Fig. 2. The geometric phase matching for (a) a single pump field and (b) two pump fields. The angle between pump fields kP1 and kP2 is θ0. The angles between the probe field kp and the pump fields kP1 and kP2 are θ1 and θ2, respectively. Here θ3 is the angle between the probe and the projected pump field. The frequency degenerate twin beams are generated with an exchange of the roles of the pump and probe beams.

In this section, we firstly theoretically describe the frequency non-degenerate and degenerate squeezed light based on a non-collinear 4WM process. As shown in Figs. 3(a) and 3(b), the double-Λ four-level processes of frequency non-degenerate and degenerate are the same except that the magnitudes of the detunings are different. Therefore, the frequency-degenerate and non-degenerate squeezed lights based on non-collinear 4WM can be described by a same set of equations.

Fig. 3. (a) Double-Λ scheme with two pump fields P1 and P2 for (a) the frequency non-degenerate and (b) frequency degenerate squeezed lights. States |3⟩ and |4⟩ are orthogonal linear combinations of magnetic states of the excited hyperfine levels. ΩP1 and ΩP2 are the Rabi frequencies, Δ1 (=ωP1ω31) and Δ2 (=ωP2 + δω42) are the detunings, and δ is the two-photon detuning.

We assume that the two pump fields EP1 and EP2 couple the transitions |1⟩ → |3⟩ and |2⟩ → |4⟩, respectively. The probe field couples the transition |2⟩ → |3⟩, and the conjugate field couples the transition |1⟩ → |4⟩. The transitions |1⟩ → |2⟩ and |3⟩ → |4⟩ are not dipole allowed transitions. In the dipole and rotating wave approximations, the Hamiltonian of the atoms combined with the Hamiltonian of the light-atom interaction is given by

where

Here ωn1 = ωnω1 (n = 2,3,4), σnm = |n ⟩ ⟨m| (n,m = 1,2,3,4), and are the Rabi frequencies, gn (n = p,c) are the atom-field coupling constants, and and are the slowly varying envelope operators of the probe and conjugate field.

The equations for the atomic operators σnm (n,m = 1,2,3,4) in the Heisenberg picture are given in the Appendix. Using the atomic operators to evaluate the linear and nonlinear components of the polarization at ωp and ωc, the polarization of the atomic medium at a particular frequency is given by and , where N is the number density of the atomic medium. The polarizations of the medium at frequency ωn (n = p,c) are given by

Here the two coefficients χpp and χcc describe the effective linear polarization processes for the probe and conjugate fields, respectively, and unlike the usual linear coefficients, they depend nonlinearly on the pump field. The other two coefficients χpc and χcp are responsible for the 4WM process. A detailed calculation is given in the Appendix.

Under the condition of the slowly varying amplitude approximation, considering nearly co-propagating beams along the z axis, these field equations in the co-moving frame are written as

where Δkz is the projection of the geometric phase mismatch Δk = kP1 + kP2kpkc on the z axis. The solutions to the propagation equations (6) and (7) with a medium of length L are given by

where

with

The number operators of the probe beam and conjugate beam are defined as and , respectively. From the above result, we define the gain of the probe beam in the 4WM process as

where the initial condition is and . The 4WM generates a correlated probe and conjugate beams, and the relative intensity fluctuations are reduced for the amplification process. After the 4WM, the relative intensity fluctuation is given by

Hence the beams have been amplified without increasing the relative intensity noise, and they are relative intensity squeezed. The standard quantum limit (SQL) is a differential measurement equal to the total optical power, that is,

The noise figure of the process (or “degree of squeezing”) is the ratio of the measured noise to the corresponding shot-noise level for equal optical power. The typical noise figure is quoted as the noise in decibels relative to the SQL.

3. Phase matching angles

In this section, we describe how to obtain the optimal angles θ1 and θ2 between the probe field and two non-collinear pump fields when the angle θ0 between pump fields P1 and P2 is given.

As shown in Fig. 4(a), when two pump fields EP1 and EP2 are incident at an angle, the total projection of the wavevector of the pump fields onto the z-axis is 2kPZ and becomes smaller; i.e., 2 |kPZ| < |kP1| + |kP2|. The geometric phase matching condition (GPMC) is given by

where θ3 is the angle between the probe and the projected pump field. In fact, if the 4WM is efficient, the GPMC of Eq. (15) may not be satisfied, but the effective phase matching condition (EPMC) must be met:

where the refractive index , and . For the case of two non-collinear pump fields, equation (16) can be written as

Fig. 4. (a) The wavevectors kP1 and kP2 are projected onto the axis as 2kPZ. The angles between the kPZ and the pump fields kP1 and kP2 are θz1 and θz2, respectively. (b) The configuration in which the geometric phase matching condition (GPMC) is fulfilled. The configuration in which the effective phase matching condition (EPMC) is fulfilled for an effective index of refraction of the probe (c) np > 1 and (d) np < 1, with necessary geometric mismatches of (c) Δkz > 0 and (d) Δ kz < 0, respectively.

Using the conservation of energy condition ωp + ωc = ωP1 + ωP2 = 2ω0, where ω0 is the frequency of the pump field, and considering θz1 = θz2 = θ0/2, equation (17) is written as

where nc ≃ 1 due to the conjugate field with a large detuning. For a given angle θ0, when np = 1, the EPMC of Eq. (18) imposes θ3 = θ0/2. Under this condition, the GPMC Δkz = 0 is also satisfied, which is the phase matching condition in free space, where the beams are required rigorously co-propagating as shown in Fig. 4(b).

When np > 1, the EPMC of Eq. (18) is established to require θ3 > θ0/2, which means that the GPMC of Eq. (15) cannot be satisfied and will occur Δkz > 0, as shown in Fig. 4(c). Considering θ1 = θ2 = θ and using the law of cosines we obtain the angle requirement between the probe field and the pump fields

If np < 1, similarly, the EPMC of Eq. (18) requires θ3 < θ0/2 and imposes Δkz < 0, as shown in Fig. 4(d). In addition, the generated probe and conjugate beams have separate directions, which requires that the angle θ3 > 0. Furthermore, using the minimum value of the refractive index min(np) according to Eq. (18), we obtain

where (θ3)min = cos−1 [2ω0 cos (θ0/2)/(min(np) ωp + ωc)] and correspondingly,

where θmin = cos−1 [ω0 (1 + cosθ0)/(min (np) ωp + ωc)]. For the case of collinear pump fields, the angle θ0 here is equivalent to 0, where the condition θ3 < θ0/2 cannot be satisfied because θ3 cannot be less than 0. That is, when np < 1, the EPMC cannot be satisfied for 4WM driven by a single pump field. Compared to the collinear pump fields case, this is a new region.

For degenerate case, the form of Eqs. (15)–(17) is the same except for the magnitude of the wave vectors. With two strong beams with the frequency of the probe and conjugate beams, and along the direction of them, and a week beam having the frequency and direction of the previous pump, we can generate the frequency degenerate and spatial nondegenerate twin beams by changing the detunings of Δ1 and Δ2. Similar to the non-degenerate case, if np + nc = 2, then GPMC Δkz = 0 will be satisfied. If np + nc > 2 or np + nc < 2, then the corresponding GPMC Δkz > 0 or Δkz < 0 will also be obtained.

4. Gain comparison

In this section, we numerically analyze the gain of the frequency non-degenerate squeezed light, and describe the phase matching between the np > 1 and np < 1 by an angle adjustment.

In our non-degenerate experiment,[56] the state |g,m⟩ (or state |1,2⟩) involves the hyperfine levels |5S1/2, F = 2,3⟩, where the hyperfine splitting of the ground state is ω21 = 2π × 3.035 GHz, and the excited state |e⟩ (or state |3,4⟩) is |5P1/2⟩ has an excited state decay rate of γ = 2π × 5.75 MHz. The pump field is blue-detuned approximately 1 GHz to the D1 line of Rb-85 5S1/2 → 5P1/2. The powers of the pump fields EP1 and EP2 are set to 350 mW, and their waists at the crossing point are 622 μm and 596 μm, respectively. The Rabi frequencies of ΩP1 and ΩP2 are ΩP1 ≃ 28γ and ΩP2 ≃ 30γ for an effective electric dipole d = 1.47 × 10−29 Cm.[62] The atomic number density of Rb-85 at 125 °C is approximately N ≃ 4.5 × 1018 m−3.

Figures 5(a) and 5(b) show the direct susceptibilities χpp and χcc for the probe and conjugate fields as a function of the two-photon detuning δ/γ, and we obtain that χcc is far less than χpp due to large detuning. In Fig. 5(a), when δ < 0, the real part of χpp is effectively responsible for the index of refraction of the probe for the single pump field case.[59,57] However, for our two-pump-fields input case, the phase matching condition can also be satisfied on the other side δ > 0.

Fig. 5. The direct and cross susceptibilities for the probe and conjugate fields as a function of the two-photon detuning δ/γ. The solid lines are the real parts, and the dashed lines are the imaginary parts. The excited state decay rate is γ = 2π × 5.75 MHz, and γc = 0.5γ. The hyperfine splitting of the ground state is ω21 = 2π × 3.035 GHz. The detuning of pump 1 is Δ1 = 174γ. The Rabi frequencies of ΩP1 and ΩP2 are ΩP1 = 28γ and ΩP2 = 30γ, respectively.

According to Eq. (12), we plot the probe gain Gp as a function of the two-photon detuning δ and the geometric phase mismatch Δkz in the presence of a single pump field and two pump fields, as shown in Fig. 6. One can see that the maximum gains are obtained on the side δ < 0 with Δkz > 0, for both the cases. When Δkz < 0, the probe gain Gp does not exist for the single pump field input case and only occurs for the two-pump-fields input case. The bandwidth of the probe gain is relatively large due to the slow change in the refractive index with the two-photon detuning.

Fig. 6. Theoretical output probe gain Gp as a function of the two-photon detuning δ/γ and the geometrical phase match Δkz with (a) a single pump field and (b) two pump fields. Here the excited state decay rate is γ = 2π × 5.75 MHz, the decoherence rate is γc = 0.5γ, the atom density is N = 4.5 × 1018 m−3, the length of the medium is L = 12.5 mm and the pump Rabi frequencies are ΩP1 = 28γ and ΩP2 = 30γ.

The theoretical output probe gain Gp as a function of the two-photon detuning δ/γ and the probe-pump angle θ with (a) a single pump field and (b) two pump fields is shown in Fig. 7, where we consider kP1 = kP2 = kP and θ1 = θ2 = θ. It can be seen that for a single pump field, the gain and 4WM process are on the δ < 0 side as the angle θ increases due to phase matching. For the case of two pump fields, the gain and 4WM process can be achieved on the left side (line L2) or the right side (line L0 or line L1) by choosing the angle θ between the probe field and the pump fields when the angle θ0 is given.

Fig. 7. Theoretical output probe gain Gp as a function of the two-photon detuning δ/γ and the probe-pump angle θ with (a) a single pump field and (b) two pump fields. The area intersecting the dashed line is the area selected by our experimental parameters. The angle between pump fields P1 and P2 is θ0/2 = 0.615°. The angles between the probe field and the pump fields P1 and P2 are θ1 = θ2 = θ = 0.861°.

The area intersecting the dashed line L1 in Fig. 7(b) is the area selected by our experimental parameters, where θ0/2 = 0.615°, and θ = 0.861°. According to the minimum value min(Re(χpp)) = −1.609 × 10−4 in Fig. 6, using Eq. (21), we obtain

If we only choose the angle θ based on the bandwidth, we can choose line L0 [Fig. 7(b)] because it has the largest bandwidth. However, in the experiment, the angle is finely adjusted according to the degree of squeezing, and the optimum value of the angle θ is different. If the angle θ is chosen as 1.2° of line L2 in Fig. 7(b), the 4WM process driven by two pump fields will also be observed on the δ < 0 side due to the large gain. On the other hand, when the probe-pump angle θ is 1.2°, the 4WM processes independently generated by the pump field EP1 and pump field EP2 can occur and two conjugate fields c1 and c2 can also be observed, because of their gains as shown in Fig. 7(a). Therefore, the strong conjugate field c and two weak conjugate fields c1 and c2 may all occur under the condition of θ ∼ 1.2° and δ < 0. That is a two-pump forward phase-matching geometry and every single-pump phase-matching condition is satisfied.[60] However, in this region, the absorption is also large, which will affect the degree of squeezing of the two generated beams.

Here θ0 has a fundamental effect on wavevector matching in the new 4WM process, thus opening up a region in which high-intensity-difference squeezed light can be obtained over a wide bandwidth with low loss and moderate gain. The gain curve in Fig. 7(b) shifts upward as the angle θ0 increases, because the minimum value of the angle θ is greater than the angle θ0/2.

The gain of the probe field and the squeezing as a function of the two-photon detuning δ is shown in Fig. 8, where the square represents the experimental data from Ref. [56] and the solid line is a theoretical simulation. The theoretical simulations and experimental data of the gain of the probe field are in good agreement, as shown in Fig. 8(a). In the experiment, the intensity-difference-squeezed light for the nondegenerate case with quantum noise almost 7 dB reduction was generated.[56] The squeezing degree is affected by the spatial mode mismatch, optical absorption by the atomic system, optical loss in the light path, and atomic decoherence. Figure 8(b) shows the theoretical simulations and experimental data of squeezing, where the theoretical squeezing curve is reduced by 0.56 times and the agreement is not very good because these effects are not included in our model in order to clarify the physics picture concisely. The bandwidth of the gain is relatively larger than that for the single pump field case, which is advantageous for realizing wide-bandwidth frequency-degenerate and nondegenerate intensity-squeezed light. These light fields can be widely used in quantum information and other fields.

Fig. 8. The gain of the probe field (a) and the squeezing (b) as a function of the two-photon detuning δ. The parameters are as follows: γ = 2π × 5.75 MHz, γc = 0.5γ, N = 4.5 × 1018 m−3, L = 12.5 mm, ΩP1 = 28γ, ΩP2 = 30γ, θ0/2 = 0.615°, and θ1 = θ2 = 0.861°.
5. Conclusion

We have studied that a two-mode squeezed light is generated from a 4WM process driven by two pump fields crossing a small angle, where the twin beams are generated with a new phase matching condition. Different from 4WM realized by a single pump field where the gain peak can only be achieved on the δ < 0 side, the new 4WM process is implemented from the δ < 0 side to the δ > 0 side by an angle adjustment. The refractive index of the corresponding probe field np can be converted from np > 1 to np < 1, which can also be used to convert between slow light[63] and fast light.[64] On the np < 1 side, the refractive index np changes slowly with the two-photon detuning δ over a large range, which leads to a relatively large gain bandwidth. With an exchange of the roles of the pump and probe beams, the frequency degenerate and spatial nondegenerate twin beams can be generated. This type of twin beams can be combined and interfered directly on the beam splitter, and can be applied in quantum information and quantum metrology.

Reference
[1] Giovannetti V Lloyd S Maccone L 2004 Science 306 1330
[2] Giovannetti V Lloyd S Maccone L 2011 Nat. Photon. 5 222
[3] Degen C L Reinhard F Cappellaro P 2017 Rev. Mod. Phys. 89 035002
[4] Abbott B P Abbott R Abbott T D et al. (LIGO Scientific Collaboration and Virgo Collaboration) 2016 Phys. Rev. Lett. 116 061102
[5] Abbott B P Abbott R Abbott T D et al. (LIGO Scientific Collaboration and Virgo Collaboration) 2019 Phys. Rev. 9 031040
[6] Yurke B McCall S L Klauder J R 1986 Phys. Rev. 33 4033
[7] Jing J Liu C Zhou Z Ou Z Y Zhang W 2011 Appl. Phys. Lett. 99 011110
[8] Hudelist F Kong J Liu C Jing J Ou Z Y Zhang W 2014 Nat. Commun. 5 3049
[9] Chen B Qiu C Chen S Guo J Chen L Q Ou Z Y Zhang W 2015 Phys. Rev. Lett. 115 043602
[10] Qiu C Chen S Chen L Q Chen B Guo J Ou Z Y Zhang W 2016 Optica 3 775
[11] Linnemann D Strobel H Muessel W Schulz J Lewis-Swan R J Kheruntsyan K V Oberthaler M K 2016 Phys. Rev. Lett. 117 013001
[12] Lemieux S Manceau M Sharapova P R Tikhonova O V Boyd R W Leuchs G Chekhova M V 2016 Phys. Rev. Lett. 117 183601
[13] Manceau M Leuchs G Khalili F Chekhova M 2017 Phys. Rev. Lett. 119 223604
[14] Anderson B E Gupta P Schmittberger B L Horrom T Hermann-Avigliano C Jones K M Lett P D 2017 Optica 4 752
[15] Gupta P Schmittberger B L Anderson B E Jones K M Lett P D 2018 Opt. Express 26 000391
[16] Du W Jia J Chen J F Ou Z Y Zhang W 2018 Opt. Lett. 43 1051
[17] Plick W N Dowling J P Agarwal G S 2010 New J. Phys. 12 083014
[18] Ou Z Y 2012 Phys. Rev. 85 023815
[19] Marino A M Corzo Trejo N V Lett P D 2012 Phys. Rev. 86 023844
[20] Li D Yuan C H Ou Z Y Zhang W 2014 New J. Phys. 16 073020
[21] Gabbrielli M Pezzé L Smerzi A 2015 Phys. Rev. Lett. 115 163002
[22] Chen Z D Yuan C H Ma H M Li D Chen L Q Ou Z Y Zhang W 2016 Opt. Express 24 17766
[23] Sparaciari C Olivares S Paris M G A 2016 Phys. Rev. 93 023810
[24] Li D Gard B T Gao Y Yuan C H Zhang W Lee H Dowling J P 2016 Phys. Rev. 94 063840
[25] Szigeti S S Lewis-Swan R J Haine S A 2017 Phys. Rev. Lett. 118 150401
[26] Gong Q K Hu X L Li D Yuan C H Ou Z Y Zhang W 2017 Phys. Rev. 96 033809
[27] Giese E Lemieux S Manceau M Fickler R Boyd R W 2017 Phys. Rev. 96 053863
[28] Hu X L Li D Chen L Q Zhang K Zhang W Yuan C H 2018 Phys. Rev. 98 023803
[29] Wei C P Hu X Y Yu Y F Zhang Z M 2016 Chin. Phys. 25 040601
[30] Gong Q K Li D Yuan C H Ou Z Y Zhang W P 2017 Chin. Phys. 26 094205
[31] Slusher R E Hollberg L W Yurke B Mertz J C Valley J F 1985 Phys. Rev. Lett. 55 2409
[32] Schnabel R 2017 Phys. Rep. 684 1
[33] Hemmer P R Katz D P Donoghue J Cronin-Golomb M Shariar M S Kumar P 1995 Opt. Lett. 20 982
[34] Lukin M D Hemmer P R Löffler M Scully M O 1998 Phys. Rev. Lett. 81 2675
[35] Zibrov A S Lukin M D Scully M O 1999 Phys. Rev. Lett. 83 4049
[36] Balic V Braje D A Kolchin P Yin G Y Harris S E 2005 Phys. Rev. Lett. 94 183601
[37] Kolchin P Du S Belthangady C Yin G Y Harris S E 2006 Phys. Rev. Lett. 97 113602
[38] Thompson J K Simon J Loh H Vuletic V 2006 Science 313 74
[39] van der Wal C H Eisaman M D André A Walsworth R L Phillips D F Zibrov A S Lukin M D 2003 Science 301 196
[40] McCormick C F Boyer V Arimondo E Lett P D 2007 Opt. Lett. 32 178
[41] McCormick C F Marino A M Boyer V Lett P D 2008 Phys. Rev. A. 78 043816
[42] Zhang Z Wen F Che J Zhang D Li C Zhang Y Xiao M 2015 Sci. Rep. 5 15058
[43] Zhang D Li C Zhang Z Zhang Y Zhang Y Xiao M 2017 Phys. Rev. 96 043847
[44] Liu C Jing J Zhou Z Pooser R C Hudelist F Zhou L Zhang W 2011 Opt. Lett. 36 2979
[45] Ma R Liu W Qin Z Su X Jia X Zhang J Gao J 2018 Opt. Lett. 43 1243
[46] Wu M C Schmittberger B L Brewer N R Speirs R W Jones K M Lett P D 2019 Opt. Express 27 4769
[47] Boyer V Marino A M Lett P D 2008 Phys. Rev. Lett. 100 143601
[48] Boyer V Marino A M Pooser R C Lett P D 2008 Science 321 544
[49] Qin Z Cao L Wang H Marino A M Zhang W Jing J 2014 Phys. Rev. Lett. 113 023602
[50] Qin Z Cao L Jing J 2015 Appl. Phys. Lett. 106 211104
[51] Du W Kong J Jia J Ming S Yuan C H Chen J F Ou Z Y Mitchell M W Zhang W arXiv: 2004.14266
[52] Corzo N V Marino A M Jones K M Lett P D 2012 Phys. Rev. Lett. 109 043602
[53] Embrey C S Turnbull M T Petrov P G Boyer V 2015 Phys. Rev. 5 031004
[54] Wang L H Fabre C Jing J 2017 Phys. Rev. 95 051802
[55] Cao L Qi J Du J Jing J 2017 Phys. Rev. 95 023803
[56] Jia J Du W Chen J F Yuan C H Ou Z Y Zhang W 2017 Opt. Lett. 42 4024
[57] Turnbull M T Petrov P G Embrey C S Marino A M Boyer V 2013 Phys. Rev. 88 033845
[58] Corzo N Marino A M Jones K M Lett P D 2011 Opt. Express 19 21358
[59] Jasperse M Turner L D Scholten R E 2011 Opt. Express 19 3765
[60] Knutson E M Swaim J D Wyllie S Glasser R T 2018 Phys. Rev. 98 013828
[61] Davis W V Kauranen M Nagasako E M Gehr R J Gaeta A L Boyd R W 1995 Phys. Rev. 51 4152
[62] Steck D A https://steck.us/alkalidata/
[63] Boyer V McCormick C F Arimondo E Lett P D 2007 Phys. Rev. Lett. 99 143601
[64] Glasser R T Vogl Ulrich Lett Paul D 2012 Phys. Rev. Lett. 108 173902